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Radiative nucleon-deuteron capture and two- and three-body photo disintegration of the three- 
nucleon bound state are described. The description uses the purely nucleonic charge-dependent CD- 
Bonn potential and its coupled-channel extension CD Bonn + A. The A-isobar excitation yields 
an effective three-nucleon force and effective two- and three-nucleon currents besides other A-isobar 
effects; they are mutually consistent. Exact solutions of three-particle equations are employed for the 
initial and final states of the reactions. The current has one-baryon and two-baryon contributions 
and couples nucleonic with A-isobar channels. A-isobar effects on the observables are isolated. 
Shortcomings of the theoretical description are discussed and their consequence for the calculation 
of observables is estimated. 
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I. INTRODUCTION 

Photo reactions in the three-nucleon system are de- 
scribed. The available scattering energy stays below 
pion-production threshold. The description allows for 
the excitation of a nucleon to a A isobar. The excita- 
tion of the A isobar remains virtual due to the available 
energy. The A isobar is therefore considered a stable 
particle; it yields an effective three-nucleon force and ef- 
fective exchange currents besides other A-isobar effects. 
The effective three-nucleon force simulates the two-pion 
exchange Fujita-Miyazawa force Q] and the three-pion 
ring part of the Illinois forces j3| in a reducible energy- 
dependent form. The effective exchange currents are of 
two-nucleon and three-nucleon nature. Since the effec- 
tive nucleonic forces and currents are built from the same 
two-baryon coupled-channel potential and from the cor- 
responding one-baryon and two-baryon coupled-channel 
current, they are consistent with each other. Since the 
two-baryon coupled-channel potential is based on the sin- 
gle exchanges of the standard isovector mesons pi (ir) and 
rho (p) and of the isoscalar mesons omega (u) and sigma 
(a), the same meson exchanges are contained in the ef- 
fective nucleonic forces and in the effective nucleonic cur- 
rents. E.g., besides the 7r exchange of Refs. 0, also p 
exchange is included in forces and currents. 

The exact solution of the three-particle scattering 
equations is used for the description of the initial- 
and final-state interactions; the coupled-channel formu- 
lation for nucleon-deuteron scattering is developed in 
Refs. [113.01; radiative nucleon-deuteron capture and 
electromagnetic (e.m.) two-body breakup of the three- 
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nucleon bound state are described in Ref. 0L Whereas a 
separable expansion of the Paris potential @] is used in 
those early calculations, Ref. solves the three-particle 
scattering equations exactly by Chebyshev expansion of 
the two-baryon transition matrix as interpolation tech- 
nique; that technique is found highly efficient and sys- 
tematic. In this paper, the technique of Ref. Q is also 
used for the description of photo processes in the three- 
nucleon system. In contrast to Ref. ja| the underlying 
purely nucleonic reference potential is CD Bonn 9]. Fur- 
thermore, the coupled-channel extension of CD Bonn, 
called CD Bonn + A and employed in this paper, is fitted 
in Ref. ^3] to the experimental two-nucleon data up to 
350 MeV nucleon lab energy; it is as realistic as CD Bonn. 
Thus, this paper updates our previous calculations 0] of 
trinucleon photo reactions. Compared to Ref. @, the 
description is extended to higher energies, and three- 
nucleon breakup is also included; however, energetically 
the description remains below pion-production thresh- 
old. An alternative description of e.m. processes in 
the three-nucleon system is given in Refs. [Til IT^ . \Tfy : 
Refs. fill II^ . ITsT l employ a different two-nucleon poten- 
tial, an explicit irreducible three-nucleon force and a dif- 
ferent e.m. current; nevertheless, the theoretical predic- 
tions of Refs. [ill Il3l| and of this paper will turn out 
to be qualitatively quite similar. 

Section |n] recalls our calculational procedure and es- 
pecially stresses its improvements. Section IIIII presents 
characteristic results for observables; A-isobar effects on 
those observables are isolated. Section IWl discusses the 
technical shortcomings of the given results. Section 
gives a summary and our conclusions. 



II. CALCULATIONAL PROCEDURE 

The calculational procedure, including the notation, is 
taken over from Ref. [fj. We remind the reader shortly of 



2 



that procedure in order to point out changes and to de- 
scribe the extension to three-body photo disintegration, 
not discussed in Ref. @. 



Nonrelativistic model for the electromagnetic 
and hadronic interaction of baryons 
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The e.m. current acts in a baryonic Hilbert space with 
two sectors, i.e., one sector being purely nucleonic and 
one in which one nucleon (N) is turned into a A iso- 
bar. The current operator is employed in its Fourier- 
transformed form J M (Q) and in a momentum represen- 
tation, based on the Jacobi momenta (pqK) of three 
particles in the definition of Ref. |j, i.e., 



(p'q'K'|J At (Q)|pqK> = S(K' - Q K) 

x (p'qV(Q,K + )|pq) 



(1) 



In Eq. Q Q is the three-momentum transfer by the pho- 
ton; it will take on particular values depending on the 
considered reaction; in the photo reactions of this pa- 
per it is given by the photon momentum k 7 . A total- 
momentum conserving <5-function is split off; the remain- 
ing current operator J M (Q, K + ) only acts on the internal 
momenta of the three-baryon system with a parametric 
dependence on the combination K + = K' + K of total 
momenta. Since all meson degrees of freedom are frozen, 
the operator has one-baryon and many-baryon pieces. 
Besides the standard nucleonic-current part there are ad- 
ditional parts involving the A isobar which then make 
effective two- and three-nucleon contributions to the ex- 
change current, the contributions being consistent with 
each other. We take one-baryon and two-baryon contri- 
butions into account, shown in Figs. Q]- El and described 
in detail in the respective figure captions. The explicit 
forms of the considered contributions are collected in Ap- 
pendix El The horizontal lines in the diagrams indicate 
that the meson exchanges are instantaneous. The dom- 
inant meson-exchange contributions arise from tt and p 
exchanges; note, that those are the only contributions of 
two-bar yon nature taken into account in the calculations 
of Refs. [ill ll2L ll^ | . In our calculations also the nondiag- 
onal pirj and uirj contributions are taken into account 
for the currents of Figs. and El The current of Fig. 
couples purely nucleonic states with states containing one 
A isobar. In contrast to Ref. Q, the contributions be- 
tween A-isobar states of one- and two-baryon nature are 
kept as shown in Fig. though the corresponding two- 
baryon contributions will turn out to be quantitatively 
entirely irrelevant; we therefore take only the diagonal 
7r contribution into account. The current is derive d by 
the extended S-matrix method of Refs. |l 4 ll "J. lift Il7|: 
however, it satisfies current conservation only approxi- 
mately with the corresponding it and p exchanges in the 
employed two-baryon interaction Hj of CD Bonn and 
CD Bonn + A. The spatial current is systematically ex- 
panded up to first order in fc/mjv, k being a characteristic 



FIG. 1: One- and two-baryon processes contained in the 
used e.m. current. In this figure only the purely nucleonic 
processes are depicted; the nucleon is indicated by the thin 
solid line, the photon by the wavy line, and the instantaneous 
meson exchange by the dashed line. In nonrelativistic order 
the one-nucleon process contributes to the charge density and 
to the spatial current, the two-nucleon processes only to the 
spatial current. The diagonal isovector tt and p exchanges are 
taken into account in the two-nucleon processes as well as the 
nondiagonal pn-y and umy contributions. 
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FIG. 2: One- and two-baryon processes contained in the used 
e.m. current. In this figure processes are depicted in which 
one nucleon is turned into a A isobar, indicated by a thick line. 
The hermitian adjoint processes are taken into account, but 
are not diagrammatically shown. In nonrelativistic order the 
one-baryon and two-baryon processes contribute only to the 
spatial current. In the one-baryon current only the magnetic 
dipole transition is kept. The diagonal isovector tt and p 
exchanges are taken into account in the two-baryon processes 
as well as the nondiagonal pn-y and ujn-y contributions. 



baryon momentum and mpf the nucleonic rest mass. The 
charge density is used in zeroth order in k/rriN for the 
standard calculations of Sec. IIIH even photo reactions 
require the charge density operator, i.e., for the Siegert 
form of the current. 

In the perturbative spirit for the evolution of photo 
processes, the e.m. interaction Hf m - acts only once, 
whereas the hadronic interaction Hj has exactly to be 
taken into account up to all orders. We use hadronic 
channel states, seen in the initial and final states \iPi) 
and |/P/) of the photo reactions with total momenta Pi 
and P/ in the form 



|4bK) 
|*a(qKK) 
|$o(pq)^oK) 



|S)|K), 

|0 Q (qH)|K), 
|to(pqH)|K) 



(2a) 
(2b) 
(2c) 
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FIG. 3: One- and two-baryon processes contained in the used e.m. current. In this figure processes are depicted which connect 
states with a A isobar. In nonrelativistic order the one-baryon process contributes to charge density and spatial current, the 
two-baryon processes only to the spatial current. Only the diagonal isovector tt exchange is taken into account in the two-baryon 
processes. 



with the energies 



E B (K) = E B 
E a {qK) 



ed 



So(pqK) = ^~ 

m N 



K 2 

6to AT ' 
3q 2 



K 2 



4ttin 
3q 2 



6to AT ' 
K 2 



4to n 6toat' 



(3a) 
(3b) 
(3c) 



to at, ed and Eb being the average rest mass of the nu- 
cleon, the deuteron and the trinucleon binding energies; 
in contrast to the notation of Ref. @, but consistent 
with our notation of hadronic reactions 0, 0] , the rest 
mass of three nucleons is removed from the energies of 
Eqs. Q). The internal trinucleon bound state is \B), 
which is normalized to 1. The product nucleon-deuteron 
and breakup channel states in the three-nucleon cm. 
frame are \4> a (q)v a } and |</>o(pq)^o) in the notation of 
Ref. 0|, v a and vq denoting all discrete quantum num- 
bers. In both cases the cm. motion |K) is explicitly 
added to the internal motion; in the three-nucleon chan- 
nel with a photon, the total momentum P is different 
from the total momentum K of the three nucleons, bound 
in the trinucleon bound state \B); in the channels with- 
out photon P = K. 

The matrix elements of the e.m. interaction require 
fully correlated hadronic states, i.e., 



|$ B K) 
|*W(qHK) 



±iOG(E B (K) ±iO)\B)\K) 
±iOG(E a (qK)±iO) 

x -i=(l + P)|$ Q (qHK), 



l*o (pq)^oK) =± i OG(P (pqK)±iO) 

x i=(l + P)|$ (pq)^K) 

with the full resolvent 

G(Z) = (Z-H -H i y 1 ; 



(4a) 
(4b) 

(4c) 
(5) 



the free Hamiltonian H contains the motion of the cen- 
ter of mass, but the rest mass of three nucleons is taken 



out, consistent with Eqs. (pj; the permutation operator 
P symmetrizes the product states; the individual kinetic 
energy operators are of nonrelativistic form; they yield 
the eigenvalues of Eqs. 10. The hadronic states l|4hl) and 
ijlcjl are normalized to S functions without additional nor- 
malization factors. Since the hadronic interaction Hamil- 
tonian Hi acts on relative coordinates only, the full resol- 
vent reproduces the bound state \B) and correlates the 
scattering states only in their internal parts, i.e., 



|*W(q) l , a K) = |^ ± )(qK)|K), 



Itf^pqHK) 



|4 ±3 (pq>o>|K>. 



(6a) 
(6b) 



S matrix for three-body photo disintegration of 
the trinucleon bound state 



The ^-matrix and the spin-averaged and spin- 
dependent cross sections for radiative nucleon-deuteron 
capture and for two-body photo disintegration of the trin- 
ucleon bound state are given in Ref. We add now 
the corresponding quantities for three-body photo disin- 
tegration. The kinematics of all considered photo pro- 
cesses is shown in Fig. 31 The figure also defines the 
employed notation for the individual particle momenta 
of the trinucleon bound state, deuteron, a nucleon, the 



fci 



k. 






FIG. 4: Schematic description of all considered three-nucleon 
photo reactions. The lines for the two-baryon and three- 
baryon particles are drawn in a special form to indicate their 
compositeness. 
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three-nucleons of break-up and the photon; k B , kd, kjsr, 
ki and fc 7 are on-mass-shell four-momenta. The corre- 
sponding particle energies are the zero components of 
those momenta, i.e., k B c, k^c, k° N c, k®c and k®c; they 
are relativistic ones with the complete rest masses in con- 



trast to those of the nonrelativistic model calculation of 
baryonic states in Eqs. JBJ- 

We give various alternative forms for the S'-matrix el- 
ements: 



(fP f \S\iPi) = - i(27r£) 4 (5(fci + fc 2 + k 3 - fc 7 - feu) (sf\M\si) (2tt?i)- 15/2 [2k°c2k° B c2k%c2k%c2k$c 
{fP f \S\iPi) = - 2mS(E N (k 1 ) + E N (k 2 ) + E N (k 3 ) - k°c-E B {k B )) (fP f \iOG(E t + iQ)Hj m \iP, 
(fP f \S\Wi) = - 2m S(E N (ki) + E N (k 2 ) + E N (k 3 ) - fc°c - E B {k B )) *(ki + k 2 + k 3 - k 7 - k B ) 



-1/2 



(4tt) 1 /2^ 



<V^ ) ( P/ q / )^o/|j' 1 (k 7 ,K + )£ Al (k 7 A)|i?). 



(2 7 rft) 3 / 2 (2fc0 C ) 1 /2 



(7a) 
(7b) 

(7c) 



Equation Il7a|) introduces a covariant form, whereas 
Eqs. l|7b|l and f7cl are noncovariant quantum mechan- 
ical realizations of it. e M (k 7 A) is the polarization vec- 
tor of the real photon with helicity A. (sf\M\si) is the 
singularity-free matrix element for three-nucleon photo 
disintegration, from which the differential cross section 



dai^f = \ {sf\M\si 



,2dLips(A; 7 + k B , fc 1; k 2 , k 3 ) 
4c 2 fc r • k^ 



(8) 



is obtained. Its dependence on the helicity A of the pho- 
ton and on the spin projection M B of the trinucleon 
bound state in the initial channel, collectively described 
by Si, and on the spin projections m s , of nucleons in the 
final channel, collectively described by Sf, are explicitly 
indicated. (sf\M \si) is Lorentz-invariant in a relativistic 
description and can therefore be calculated in any frame. 
However, in our model it is calculated in the framework 
of nonrelativistic quantum mechanics and therefore loses 
the property of being a Lorentz scalar; equating Eqs. ifTaj) 
and ijTcjl (sf\M\si) is defined by 



(s f \M\ Si 



(2nh) 3 [2k° B c 2fc?c 2k%c 2fcgc] 1/2 

x (4- ) (p / q / )^ 0/ |^(k 7 ,K + )e, 1 (k 7 A)|i?). 

(9) 



We calculate that matrix element in the center of mass 
(cm.) system of the final hadronic state using the follow- 
ing computational strategy. The strategy is nonunique, 
since the model calculations, due to dynamic limitations, 
miss the trinucleon binding energy; the necessary correc- 
tion for that miss has arbitrary features. In contrast, 
the S matrix of Eq. itTajl is based on proper relativistic 
kinematics with experimental rest masses. 

1. The experimental photon momentum k 7 in the 
lab frame with k^ = determines the total momentum 
Pf and energy i? (p/q/K^) of the final three-nucleon 
system in the lab frame, i.e., P/ = Ky = k 7 and 



the experimental trinucleon binding energy. The result- 
ing energy i?o(p/q/K/) of the final state is the true ex- 
perimental one. Thus, the experimental two-body and 
three-body breakup thresholds are exactly reproduced. 
In nonrelativistic approximation for baryon kinematics, 
the internal three-nucleon kinetic energy part of the final 
state is p^/niN + 3q 2 r/4mAr = i?o(p/q/K/) — K 2 /6mjv. 

2. The matrix element (sf\M\si) is calculated in the 
cm. system as on- energy- shell element under nonrel- 
ativistic model assumptions. Under those assumptions 
the internal energy of the initial state is |k 7cm |c + E B + 
k 7cm /6mjv = i?o(p/q/K/) - Kj/6m N , k 7cm being 
the photon momentum in the cm. system, in which 
the trinucleon bound state is moving with momentum 
k B = — k 7cm ; thus, K + = — k 7cm . Taking the computed 
trinucleon model binding energy E B and the average nu- 
cleon mass m^v, i.e., mjyc 2 = 938.919 MeV, the magni- 
tude of the photon momentum |k 7cm | to be used for the 
current matrix element results. Since the model binding 
energy E B is not the experimental one, neither for 3 He 
nor for 3 H, and since the cm. contribution to total three- 
nucleon energies is assumed to be nonrelativistic with 
mass 3m n and to separate from its internal part, that 
photon momentum k 7 cm does not have the experimental 
value. In the matrix element (sf\M\si), to be calculated 
according to Eq. ©, the zero-momentum components 
kf and kP B are nonrelativistic model quantities and dif- 
fer from the baryonic energies just by rest masses, i.e., 
kfc — E N (ki) + m N c 2 and k%c = E B {k B ) + 3m N c 2 . 

In contrast to the matrix element (sf\M\si) which car- 
ries the dynamics, the kinematical factors in Eq. i.e., 
the Lorentz-invariant phase-space element 



dLips(fc 7 + fc_B,fci, k 2 , k 3 ) 



= (2itKf S(h + k 2 + k 3 ~k 1 -k B ) 
d 3 k x d 3 k 2 d 3 k 3 



(10) 



#o(p/q/K/) = E B 



|k 7 |c. 



This step is done using 



(2nh) 9 2k\c 2k\c 2k%c 
and the factor 4c 2 k B ■ k 7 , which contains the incoming 
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flux, the target density and projectile and target normal- 
ization factors can be calculated relativistically. 

The momenta in the initial and final states are con- 
strained by energy and momentum conservation. E.g., if 
the momentum ki and the direction k2 were measured, 
all three nucleon momenta are determined in the final 
state, although not always uniquely. In practice, the two 
nucleon scattering angles with respect to the beam direc- 
tion (di,(pi) and (0 2 ,(p 2 ), usually notationally shortened 
to (0i, #2 1 ^ — fi), and their kinetic energies without rest 
masses E\ — En (hi) and E<i = E^^k-i) are measured. 
Those energies are related by momentum and energy con- 
servation and therefore lie on a fixed kinematical curve. 
The observables are therefore given as function of the 
arclength S along that curve, i.e., 



dS 



(11) 



with dS — y/ dE\ + dE 2 and E% being considered a func- 
tion of E\ or vice versa depending on numerical conve- 
nience; the arclength is always taken counterclockwise 
along the kinematical curve. No confusion between the 
arclength S and the iS-matrix of Eq. @ should arise. 
The normalization of the arc length value zero is chosen 
as max{i?i|i?2 = 0}. 

The lab cross section therefore takes the compact form 



d 5 <ji-^f = 1 (s f | M | S-) | fps dS d 2 ki d 2 k 2 



(12a) 



with the abbreviation fps for a phase-space factor; in the 
lab frame fps is 



fps 



fps 



(2TTH)- 5 

6{ki + 



d 3 k 3 k\ dk 2 



dS 



fc 7 - k B ) 



(12b) 
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k 2 k 2 
K 1 K 2 



32c 7 k°m B 
x {k?[|k 2 |(fc° + fc 3 °) - k° 2 k 2 ■ (k 7 -k0] 2 

+ k2[|ki|(fc0 + ^)-fciki-(k 7 ~k 2 )] 2 } 



-1/2 

(12c) 



The cross section 112all is still spin-dependent. The spin- 
averaged fivefold differential cross section is 



d 5 a 



d 5 (T. 



dSdsii dn 2 4 j£ x msi ^ ms3 dsd 2 ki d 2 k 2 



(13) 



Spin observables are defined as in Refs. 0,0]. The ex- 
perimental setup determines the isospin character of the 
two detected nucleons 1 and 2. 

The calculational strategy of Eqs. © - ltl2ll is in the 
spirit of Ref . Q ; it chooses the kinematics differently for 



the dynamic matrix element (sf\M\si) on one side and 
for the phase space dLips(fc 7 +&B, fci, k 2l k 3 ) and the factor 
4c 2 &b 'ky on the other side. That strategy can be carried 
out with ease for the observables of exclusive processes. 
However, when total cross sections in hadronic and e.m. 
reactions or inelastic structure functions in electron scat- 
tering are calculated as described in Appendix El for the 
total photo cross section, the energy conserving 5 func- 
tion is rewritten as imaginary part of the full resolvent 
and has to be made consistent with the employed non- 
relativistic dynamics. Thus, as described in AnnendixIBI 
the split calculational strategy, developed in Ref. Q and 
so far here, cannot be carried through for total cross 
sections and inelastic structure functions; furthermore, 
as discussed in Ref. Q for the hadronic reactions, such 
a split calculational strategy would also be inconsistent 
with the fit of the underlying baryonic potentials. 

We shall therefore use nonrelativistic kinematics in the 
framework of quantum mechanics throughout. The cor- 
responding expressions, derived from quantum mechan- 
ics directly, can also be obtained formally from Eqs. © 
- iflfljl based on quantum field theory by replacing the 
hadron energy factors 2fc"c by their rest masses 2rrijC 2 
and using nonrelativistic energies for the energy conserv- 
ing (^-functions and for the definition of the kinematic 
locus. The lab cross section is constructed from the fol- 
lowing building blocks, i.e., the matrix element 



(s f \M\. 



,2x3ll/2 



(2irh) 3 [2m B c 2 (2m N c 2 ) 

X (P/qj>0/ \f (k 7 cm, -k 7 cm) 

x e (tl (k 7cm A)| J B) 



(14a) 



and the phase space factor fps of Eqs. l(T2|l which takes 
the following changed form 



4c 2 k^rriE 



dS 



5(E N (ki) + E N (k 2 ) + E N (k 3 ) - k° y c - E B ) 



[2m N c 2 f 
x <5(ki + k 2 + k 3 - k 7 ), 



(14b) 



fps = 



(2ttK)- 



32c 8 fc°77l NlTlB 



k 2 k 2 {k 2 [2|k 2 |-k 2 -(k 7 -k 1 )]' 



-k 2 [2|k 1 |-ki.(k 7 -k 2 )] 2 } 



-1/2 



(14c) 



Section IIVI will discuss the differences between the 
present fully nonrelativistic calculational scheme of cross 
sections and that of Eqs. and 11211 with some rela- 
tivists features. 



III. 



RESULTS 



We present results for spin-averaged and spin- 
dependent observables of nucleon-deuteron radiative cap- 
ture and of three-nucleon photo disintegration; results of 
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two-nucleon photo disintegration are transformed to cor- 
responding ones of radiative capture. The results are 
based on calculations derived from the purely nucleonic 
CD-Bonn potential @ and its coupled-channel exten- 
sion [lfT |. which allows for single A-isobar excitation in 
isospin-triplet partial waves. The A isobar is considered 
to be a stable particle of spin and isospin 3/2 with a 
rest mass m&c 2 of 1232 MeV. In contrast to the coupled- 
channel potential constructed previously by the subtrac- 
tion technique 0] and used in the calculations of Ref. @, 
the new one of Ref. ^3] is fitted properly to data and 
accounts for two-nucleon scattering data with the same 
quality as the original CD-Bonn potential. We describe 
first the standard calculational procedure adopting the 
strategy of Sec. Ill Bl 

The hadronic interaction in purely nucleonic and in 
nucleon-A partial waves up to the total two-baryon an- 
gular momentum I = 4 is taken into account. The cal- 
culations omit the Coulomb potential between charged 
baryons. Nevertheless, the theoretical description is 
charge dependent. For reactions on 3 He the pp and np 
parts of the interaction are used, for reactions on 3 H the 
nn and np parts. Assuming charge independence, the 
trinucleon bound state and nucleon-deuteron scattering 
states are pure states with total isospin T = the three- 
nucleon scattering states have total isospin T = \ and 
T = §, but those parts are not dynamically coupled. 
Allowing for charge dependence, all three-baryon states 
have T = \ and T = § components which are dynam- 
ically coupled. For hadronic reactions that coupling is 
found to be quantitatively important in the 1 So partial 
wave ; in other partial waves the approximative treat- 
ment of charge dependence as described in Ref. [19] is 
found to be sufficient; it does not couple total isospin 
T = 5 and § channels dynamically. The same applies for 
photo reactions considered in this paper. The effect of 
charge dependence is dominated by the 1 S'o partial wave; 
it is seen in some particular kinematics of radiative cap- 
ture and of three-body photo disintegration; we do not 
discuss it in this paper. Furthermore, the calculations of 
e.m. reactions require total isospin T = | components 
of scattering states in all considered isospin-triplet two- 
baryon partial waves, since the e.m. current couples the 
T = g and T = § components strongly. 

The three-particle equations for the trinucleon bound 
state LB) and for the scattering states are solved as in 
Ref. pj; in fact, the scattering states are calculated only 
implicitly as described in Appendix [bJ The resulting 
binding energies of 3 He are -7.941 and -8.225 MeV for CD 
Bonn and CD Bonn + A, respectively. If the Coulomb in- 
teraction were taken into account, as proper for 3 He, the 
binding energies shift to -7.261 and -7.544 MeV, whereas 
the experimental value is -7.718 MeV. Nevertheless, we 
use the purely hadronic energy values and bound-state 
wave functions for consistency when calculating the cur- 
rent matrix elements, since we are unable to include the 
Coulomb interaction in the scattering states. 

Whereas the hadronic interaction is considered up to 



1 = 4, the e.m. current is allowed to act between partial 
waves up to / = 6, the higher partial waves being cre- 
ated by the geometry of antisymmetrization. The e.m. 
current is taken over from Refs. (TS . l2fT | with some nec- 
essary modifications: 1) The e.m. current is richer than 
the one used in Ref. [g]; diagonal two-baryon currents 
connecting states with A isobar are taken into account. 
2) More recent values for the e.m. couplings of the A 
isobar are used according to Refs. [2ll l22l|. 3) Meson 
coupling constants, meson masses and hadronic form fac- 
tors used in meson-exchange currents (MEC) are chosen 
consistently with the employed hadronic interactions CD 
Bonn and CD Bonn + A; they are listed in Refs. |flLHcj. 
The employed contributions to the e.m. current are col- 
lected in Appendix El The current is expanded in elec- 
tric and magnetic multipoles as described in Refs. 
The technique for calculating multipole matrix elements 
is developed in Ref. [ID]; a special stability problem @ 
arising in the calculation requires some modifications of 
that technique as described in Ref. [^j. The magnetic 
multipoles are calculated from the one- and two-baryon 
parts of the spatial current. The electric multipoles use 
the Siegert form of the current without long-wavelength 
approximation; assuming current conservation, the dom- 
inant parts of the one-baryon convection current and of 
the diagonal n- and p-exchange current are taken into 
account implicitly in the Siegert part of the electric mul- 
tipoles by the Coulomb multipoles of the charge density; 
the remaining non-Siegert part of the electric multipoles 
not accounted for by the charge density is calculated us- 
ing explicit one- and two-baryon spatial currents. The 
charge density contributing to the Siegert term has di- 
agonal single-nucleon and single-A isobar contributions 
only; the nucleon-A transition contribution as well as 
two-baryon contributions are of relativistic order and are 
therefore omitted in the charge-density operator when 
calculating Coulomb multipoles. 

The number of considered current multipoles is limited 
by the maximal total three-baryon angular momentum 
Jmax = -y , taken into account for the hadronic scattering 
states. The results for the considered photo reactions 
up to pion-production threshold appear fully converged 
with respect to higher two-baryon angular momenta /, 
with respect to A-isobar coupling and with respect to 
higher three-baryon angular momenta J on the scale of 
accuracy which present-day experimental data require. 

That is the standard calculational procedure. Sec- 
tion IIVI describes the shortcomings of that standard de- 
scription. In the rest of this section we focus on A-isobar 
effects in sample observables. 



A. Nucleon-deuteron radiative capture 

Figures and present results for spin-averaged and 
spin-dependent observables of radiative nucleon-deuteron 
capture at 100 and 150 MeV nucleon lab energy; a result 
for 200 MeV will be shown later in Sec. II V C 31 Results 
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FIG. 5: Differential cross section and analyzing powers of proton-deuteron radiative capture at 100 MeV nucleon lab energy as 
function of the cm. nucleon-photon scattering angle. Results of the coupled-channel potential with A-isobar excitation (solid 
curves) are compared with reference results of the purely nucleonic CD-Bonn potential (dashed curves). The experimental data 
are from Ref. [24] (circles) and from Ref. [2^| (crosses). 
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FIG. 6: Differential cross section and nucleon analyzing power of proton-deuteron radiative capture at 150 MeV nucleon 
lab energy as function of the cm. nucleon-photon scattering angle. Results of the coupled-channel potential with A-isobar 
excitation (solid curves) are compared with reference results of the purely nucleonic CD-Bonn potential (dashed curves). The 
experimental data are from Ref. [2fil ]. 



for the time-reversed two-body photo disintegration of 
the trinucleon bound state are not shown separately. The 
energies are well above those of Ref. Q, but remain be- 
low pion-production threshold. Control calculations at 
lower energies indicate that the results of Ref. Q do not 
get any essential physics change, though the hadronic in- 



teraction and the e.m. current are improved compared 
with Ref. 0. 

There are noticeable A-isobar effects on the considered 
observables, especially on the nucleon analyzing power 
A y (N); it is described rather well with the inclusion of 
the A isobar. Reference |24| presents new experimental 
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FIG. 7: Total 3 H three-nucleon photo disintegration cross 
section as function of the photon lab energy Ej. Results of 
the coupled-channel potential with A-isobar excitation (solid 
curve) are compared with reference results of the purely nu- 
cleonic CD-Bonn potential (dashed curve). The experimental 
data are from Ref. [2tH |. 
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FIG. 8: The fourfold differential cross section of the 
3 He(7,pn)p reaction at 85 MeV photon lab energy as function 
of the pn opening angle at 9 P = 81°. Results of the coupled- 
channel potential with A-isobar excitation (solid curve) are 
compared with reference results of the purely nucleonic CD- 
Bonn potential (dashed curve). The experimental data are 
from Ref. H3- 



data for the differential cross section and deuteron vector 
analyzing power at 200 MeV deuteron lab energy, corre- 
sponding to 100 MeV nucleon lab e nerg y. There is a dis- 
crepancy between new |24| and old |25l | differential cross 
section data in the maximum region; the new data are in 
good agreement with our results including the A isobar. 
However, there is a clear disagreement between theoret- 
ical predictions and experimental data at small scatter- 
ing angles getting more pronounced at higher energies; 
one possible reason for that discrepancy is discussed in 
Sec. lIV (T3l There is also a modest beneficial A-isobar ef- 
fect on the deuteron vector analyzing power A y (d) . The 
theoretical prediction for one deuteron tensor analyzing 
power, i.e., A xx , is also given in Fig.0 our motivation for 
showing A xx is the fact that an experiment determining 
deuteron tensor analyzing powers is in progress 

Our results are qualitatively consistent with those of 
Refs. EJHII3- 



B. Three-body photo disintegration of 
three-nucleon bound state 



putational strategy as discussed in Sec. lIV B"T1 Figures|H| 
-Elshow semi-exclusive fourfold differential cross sections 
of 3 He photo disintegration at higher energies; they are 
obtained from the fivefold differential cross section 11 ■'•ill 
by integrating over the kinematical curve S. Again, the 
A-isobar effect for those particular observables appears 
rather small, smaller than the experimental error bars. 
There is also disagreement between theoretical predic- 
tions and experimental data in some kinematical regimes 
which in part may be caused by experimental conditions, 
e.g., by finite geometry, not taken into account in our 
calculations. 

Finally, in Fig. EH we show fully exclusive sam- 
ple fivefold differential cross sections of three-nucleon 
photo disintegration at 120 MeV photon lab energy for 
two kinematical configurations which were shown semi- 
exclusively in Fig. even at that higher energy the A- 
isobar effect is rather mild. 



Experimental data for three-nucleon breakup are much 
scarcer than for two-body photo disintegration. To the 
best of our knowledge, there are no fully exclusive experi- 
mental data in the considered energy regime; we therefore 
show in Figs.0-|H|our predictions for inclusive and semi- 
exclusive observables and compare them with existing ex- 
perimental data. Figured shows our results for the total 
3 H three-nucleon photo disintegration cross section in the 
low energy region- there is no significant A-isobar effect. 
In contrast, Ref. [lj| sees a larger three-nucleon force ef- 
fect for this observable; this discrepancy is partly due to 
a larger three-nucleon force effect on trinucleon binding 
and subsequent scaling and partly due to a different com- 



IV. SHORTCOMINGS OF THE DESCRIPTION 

The present description of photo reactions is with re- 
spect to the dynamic input, i.e., with respect to the 
hadronic interaction and to the e.m. current, and with 
respect to the scope of applications a substantial im- 
provement compared with Ref. [f|. But it is still not 
a unique and in itself consistent description. We are un- 
able to repair the existing deficiencies. However, this sec- 
tion points those shortcomings out and tries at least to 
estimate their size. We identify three different problem 
areas. 
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FIG. 9: The fourfold differential cross section of the 
3 He(7,pp)n reaction as function of the photon lab energy E y 
in various kinematical configurations: (81.0°, 81.3°, 180.0°) 
(top), (92.2°, 91.4°, 180.0°) (middle), and average of 
(81.5°, 90.8°, 180.0°) and (91.7°, 80.9°, 180.0°) (bottom). Re- 
sults of the coupled-channel potential with A-isobar excita- 
tion (solid curve) are compared with reference results of the 
purely nucleonic CD-Bonn potential (dashed curve). The ex- 
perimental data are from Ref. 28]. 



A. Shortcomings of the theoretical form of the 
cross section 

Our standard strategy uses the nonrelativistic 
form (1 1-41) for cross sections; this choice appears to be con- 
sistent with the underlying two-baryon dynamics, though 
inconsistent with the experimental relativistic kinemat- 
ics. We therefore compare results obtained from Eqs. ifl^jl 
with corresponding ones obtained from the relativistic 
form of the cross section ll"12ll which uses relativistic ki- 
netic energies for the Lorentz-invariant phase space ele- 
ment iflfljl and the kinematic locus lfTT|l combined with 
the nonrelativistic matrix element 0. The comparison 
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FIG. 10: The fivefold differential cross section of three- 
nucleon photo disintegration at 120 MeV photon lab en- 
ergy as function of the arclength S along the kinematical 
curve for configuration (92.2°, 91.4°, 180.0°) on the top and 
(81.5°, 90.8°, 180.0°) on the bottom. Results of the coupled- 
channel potential with A-isobar excitation (solid curve) are 
compared with reference results of the purely nucleonic CD- 
Bonn potential (dashed curve). 



is possible for observables in fully exclusive reactions. 

The difference between those aspects of relativistic and 
nonrelativistic kinematics is minor for all considered ob- 
servables of radiative capture, i.e., less than 1%, but more 
significant, i.e., up to 10%, for three-nucleon photo disin- 
tegration as shown in Fig. El In ah considered cases, the 
relativistic and nonrelativistic kinematical curves 11 111 
are very close to each other, e.g., for the configuration 
of Fig. El the distance between them in the E\ — E 2 
plane is 0.5 MeV at most, and their total arclengths are 
140.2 MeV and 142.0 MeV, respectively. For the compar- 
ison the nonrelativistic results in Fig. ^] are scaled down 
to the relativistic arclength by the factor 140.2/142.0; 
the shown change, however, is due to the difference in 
the phase space factors fps of Eqs. ill2bj) and l|14bjl . 

We emphasize: The effect indicated in this subsection 
does not represent the true difference between nonrela- 
tivistic quantum mechanical and fully relativistic quan- 
tum field theoretical results, but it may indicate the order 
of magnitude of the shortcomings of nonrelativistic calcu- 
lations. In the light of the accuracy of present day data, 
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FIG. 11: Differential cross section of three-nucleon photo 
disintegration at 120 MeV photon lab energy as function of 
the arclength S along the kinematical curve for configuration 
(91.7°, 80.9°, 180.0°). Results of the coupled-channel poten- 
tial with A-isobar excitation based on nonrelativistic phase 
space (solid curve) according to Eqs. 114(1 are compared with 
results based on relativistic phase space (dashed curve) ac- 
cording to Eqs. 11211 . 



this shortcoming of the theoretical description is rather 
inconsequential. 
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FIG. 12: Differential cross section of proton-deuteron ra- 
diative capture at 19.8 MeV deuteron lab energy as func- 
tion of the cm. nucleon-photon scattering angle. Results 
of the coupled-channel potential with A-isobar excitation de- 
rived from the standard approach (solid curve) are compared 
with results of option three which uses experimental energies 
for both initial and final states, but the matrix element Q14;il 
is off-shell (dashed-dotted curve). The results of option two 
are rather close to the solid curve. In order to appreciate the 
effect of the nonunique choice of kinematics in relation to the 
size of the A-isobar effect, results of a standard calculation 
with the purely nucleonic reference potential are also given as 
dashed curve. The experimental data are from Ref. |2fll |. 



B. Shortcomings of the dynamics 

1. Nonunique choice of kinematics 

Our computational strategy in choosing the kinematics 
for the matrix element (sf\M\si) is described in Sec. Ill Bl 
(sf\M\si) is calculated in the cm. system. We opt to 
let the experimental beam energy determine the energy 
of hadronic nucleon-deuteron state in radiative capture 
and the energy of the hadronic two-body and three-body 
final states in photo disintegration exactly. Since the trin- 
ucleon model binding energy is not the experimental one 
and the kinematics is nonrelativistic for baryons when 
calculating (sf\M\si), the energy of the photon does not 
have the experimental value when assuming energy con- 
servation. At very low energies the deviation can get as 
large as 10%, whereas at higher energies considered in 
this paper it remains around 1 - 2%. In contrast, in a 
second option we could let the experimental beam en- 
ergy determine the cm. photon energy exactly; then 
the energies of the hadronic nucleon-deuteron and three- 
nucleon states are not experimental ones. A third option 
may use experimental energies for both initial and final 
states, but then the matrix element determining physical 
amplitudes is slig htly off-shell; this is the computational 
strategy of Refs. [I1E1- 

The difference in results be- 
tween those three choices is minor at higher energies, 
i.e., above 100 MeV nucleon lab energy, for all consid- 
ered observables in all considered kinematical regimes. 
However, there are differences up to 10% for observables 



at low energies. There, the observed A-isobar effect de- 
pends strongly on the choice of computational strategy. 
An example is shown in Fig. 1121 



2. Omission of Coulomb interaction between protons 

We are unable to include the Coulomb interaction in 
the three-nucleon scattering states. In contrast, the se- 
lected inclusion of the Coulomb interaction in the trin- 
ucleon bound state is easily possible, but this inclu- 
sion creates an additional inconsistency: Initial and final 
hadronic states become eigenstates of different Hamilto- 
nians, and, strictly speaking, the Siegert form of the cur- 
rent operator is not applicable. Nevertheless, we do such 
an inconsistent calculation which Refs. |12l [lj| chooses 
to do as standard calculation, in order to estimate the 
effect of the omitted Coulomb interaction at least par- 
tially. The inclusion of the Coulomb interaction in the 
trinucleon bound state systematically reduces the spin- 
averaged cross sections; in contrast, spin observables ap- 
pear to be almost unaffected. A characteristic result 
is shown in Fig. Even at higher energies the ob- 
served Coulomb effect may be of the same order of mag- 
nitude as the full A-isobar effect; however, it is not clear 
whether the indicated effect represents a true Coulomb 
effect or just the inconsistency between bound and scat- 
tering states. 
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FIG. 13: Differential cross section of proton-deuteron radia- 
tive capture at 95 MeV deuteron lab energy as function of 
the cm. nucleon-photon scattering angle. Results of the 
coupled-channel potential with A-isobar excitation derived 
from the standard approach (solid curve) are compared with 
results including the Coulomb interaction in the three-nucleon 
bound state (dashed curve). The experimental data are from 
Ref. H}. 



C. Shortcomings of the e.m. current 

1. Lack of current conservation 

The potentials CD Bonn and CD Bonn + A used in 
this paper have nonlocal structures, whereas the e.m. 
current, given explicitly in Appendix El is employed in a 
local nonrelativistic form. Thus, the continuity equation 
is not fulfilled for the current. As measure for this de- 
ficiency predictions are compared based on two different 
approaches for the electric multipoles, i.e., (1) the stan- 
dard calculation with the Siegert operator accounting 
for the two-baryon currents implicitly by assumed cur- 
rent conservation and (2) the explicit use of the meson- 
exchange currents for all of the electric multipoles. The 
discrepancy between those two calculations measures the 
importance of the existing lack of current conservation; 
indeed the violation can be significant as Fig. 1141 proves. 
We believe that calculations with the Siegert form of the 
current operator repair the violation of current conserva- 
tion in part; we therefore employ the Siegert form of the 
current operator in our standard calculational strategy. 

However, at this stage it is useful to discuss the lack of 
current conservation in more detail: 

1) The a, p and ui exchanges yield a spin-orbit inter- 
action. That spin-orbit interaction makes, even in local 
approximation and even for isoscalar-meson exchanges, 
a contribution to the continuity equation. The corre- 
sponding spin-orbit contribution to the exchange current 
is derived in local form pj^ . There are also additional 
contributions to the p-exchange current not taken 
into account in our standard calculations. All that con- 
tributions are implicitly contained in the Siegert-part of 
the electric multipoles. In the tentative calculations de- 
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FIG. 14: Differential cross section of proton-deuteron radia- 
tive capture at 190 MeV nucleon lab energy as function of 
the cm. nucleon-photon scattering angle. Results of the 
coupled-channel potential with A-isobar excitation derived 
from the Siegert approach for electric multipoles (solid curve) 
are compared with results based on the explicit use of meson- 
exchange currents (dashed curve). In order to appreciate the 
size of the two-baryon current contribution, the results of a 
non-Siegert calculation with one-baryon currents only are also 
given as dashed-dotted curve. The experimental data are from 
Ref. [al- 



scribed in this paragraph they are used explicitly for the 
non-Siegert part of the electric multipoles and for the 
magnetic multipoles. There, that sample contributions 
yield only small corrections for observables considered in 
this paper, seen mostly in spin observables; however, even 
for them the corrections are of the order of 2% at most. 
We therefore conclude that spin-orbit contributions and 
additional p-exchange currents can quite safely be ne- 
glected in the Siegert form of the current operator when 
calculating the photo reactions of this paper. However, 
they are more important for calculations based fully on 
explicit exchange currents and they therefore make a non- 
negligible contribution to the difference seen in Fig. El 

2) The employed potentials have explicit nonlocal 
structures. That explicit nonlocality makes a contribu- 
tion to the two-baryon exchange current. 

We compared the results of the following models for the 
purely nucleonic potential. The models are based on n, 
p and oj exchanges and on a scalar isoscalar and a scalar 
isovector a exchange; one model uses the nonlocal struc- 
tures of CD-Bonn type and the other their local approxi- 
mations. Both models are tuned to deuteron binding and 
to 1 Sq and 3 Si — 3 Z?i phase shifts only. Though nonre- 
alistic models, both usually predict observables in qual- 
itative agreement with the realistic descriptions. When 
employing the local potential model we obtain identi- 
cal results for calculations based on the Siegert form of 
the current operator and calculations based on explicit 
exchange-current contributions to all multipoles. How- 
ever, when employing the nonlocal potential model, the 
results can differ substantially; the difference can be as 
large as seen in Fig. El f° r the coupled-channel poten- 
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tial CD Bonn + A. We conclude, the explicit nonlocal- 
ity of the employed potentials is a significant source for 
current non-conservation. Future calculations should at- 
tempt to design nonlocal exchange-current contributions 
consistent with the nonlocality of the underlying baryon- 
baryon potentials. 

3) The employed potentials have an implicit nonlocal- 
ity due to the general partial- wave dependence of the me- 
son exchanges. That partial-wave dependence is slight for 
the 7r, p and w exchanges, but substantial for the a ex- 
change. That implicit nonlocality makes a contribution 
to the two-baryon exchange current. 

The nucleonic CD Bonn as well as CD Bonn + A show 
a small partial-wave dependence in 7r and in p exchange. 
Fitting CD Bonn with partial wave independent 7r and p 
exchanges decreases the quality of the fit to data only by 
very little; % 2 /datum increases from 1.02 to 1.03. When 
comparing observables of the hadronic processes and of 
the photo reactions of this paper for both potentials, no 
distinguishable difference is found in plots. We conclude: 
The implicit nonlocality arising from the partial-wave 
dependence in 7r and p exchange of CD Bonn and CD 
Bonn + A is of no consequence for the prediction of 
observables. 

The local model used for the discussion of problem 2) 
is modified to simulate the partial-wave dependence of 
u> exchange in the nucleonic CD Bonn, it is retuned as 
under 2). The u> exchange is taken to be without hadronic 
cutoff in the l P\ partial wave as in CD Bonn; this partial- 
wave dependence violates current conservation. However, 
the observed difference between calculations based on 
the Siegert form of the current operator and calculations 
based on explicit exchange-current contributions to all 
multipoles is much smaller than that shown in Fig. El 
We conclude: The implicit nonlocality arising from the 
partial-wave dependence in the w-exchange of CD Bonn 
and CD Bonn + A is of no real consequence for the 
prediction of observables. 

With respect to the partial-wave dependence of cr ex- 
change the local model used for the discussion of problem 
2) is studied. We concentrate on the difference of a ex- 
change between isospin singlet and triplet partial waves, 
i.e., on the effective isovector nature of the a meson in- 
troduced in the model and in CD Bonn and CD Bonn 
+ A. Furthermore, even if the a exchange were truly 
an isoscalar one in purely nucleonic potential, the ex- 
plicit treatment of the A isobar in the coupled-channel 
extension introduces an isovector correction: The em- 
ployed coupled-channel potential CD Bonn + A, act- 
ing in isospin-triplet partial waves, has a weakened a 
exchange compared to the purely nucleonic CD Bonn; 
part of the intermediate range attraction simulated by 
a exchange is taken over by A-isobar excitation in the 
coupled-channel approach. Thus, the A-isobar current 
has to be supplemented by changed c-exchange current. 
Omitting the cr-meson contribution to the exchange cur- 
rent, quite significant differences, comparable to that of 
Fig- El arise for observables of the photo reactions in this 



paper between calculations based on the Siegert form of 
the current and on the full explicit exchange-current con- 
tributions. In contrast, the explicit tr-meson exchange- 
current contributions to the non-Siegert part of the elec- 
tric multipoles and to the magnetic multipoles remain 
small. We arrive to qualitatively the same results when 
including the cr-meson exchange-current for CD Bonn 
and CD Bonn + A with the cr-meson parameters of S- 
waves. We conclude: Though the partial-wave depen- 
dence of the cr-meson exchange is a significant source of 
current nonconservation, the standard calculation based 
on the Siegert form of the current for part of the elec- 
tric multipoles and on explicit exchange-current contri- 
butions to all other multipoles appears to be quite a re- 
liable calculational scheme. 



4) The employed potentials are charge dependent. The 
charge dependence of the interaction is due to the charge 
dependence of the parameters of exchanged n, p and a 
mesons and due to the charge dependence of the nu- 
cleonic masses. The isospin structure of the charge- 
dependent potential contributions is given in terms of the 
baryonic isospin projections; thus, that isospin depen- 
dence, giving rise to charge dependence, does not require 
an exchange current by itself; it only does so, if its po- 
tential forms were nonlocal. In the case of the employed 
potentials it is so indeed, but that explicit nonlocality was 
already discussed in problem 2). The diagonal tt- and p- 
exchange contributions to the exchange current should be 
built from the meson parameters of the charged mesons. 
The nondiagonal it and p exchanges are carried by the 
mesons of all charges. However, our standard calcula- 
tion uses averaged meson parameters and an averaged 
nucleon masses for all meson-exchange currents; it was 
checked that both calculational simplifications are with- 
out any consequence for the observables of this paper. 



From this lengthy, but we think necessary discussion 
of the problems 1) to 4) we conclude for the calcula- 
tions of this paper: When the Siegert form of the current 
is used for part of the electric multipoles and explicit 
exchange-current contributions to all other multipoles in 
the operator form of Appendix El the implicit nonlo- 
cality of CD Bonn and CD Bonn + A arising from the 
partial- wave dependence of the meson exchanges is with- 
out consequences for prediction. In contrast, the explicit 
nonlocality of CD Bonn and CD Bonn + A, also respon- 
sible for current nonconservation, is of serious concern; 
its consequence on the non-Siegert parts of the current 
could not be estimated yet by any of our models. Still, 
we believe that our standard calculation, based on the 
Siegert form of the current, effectively corrects the cur- 
rent nonconservation and is therefore quite reliable for 
the observables of photo reactions considered in this pa- 
per. 
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2. Lack of covariance 

If a fully covariant description of dynamics were 
available, the current matrix element (sf\M\si) were a 
Lorentz scalar and could therefore be calculated in any 
frame with identical results. However, our description of 
hadron dynamics is nonrelativistic, and the results there- 
fore are frame-dependent. We investigate that frame de- 
pendence calculating the same matrix elements in lab 
and in cm. frames, i.e., in the rest frames of the initial 
and final three- nucleon systems. The two frames differ 
by the three-nucleon total momentum and by the pho- 
ton momentum. However, we found that for the observ- 
ables considered in this paper the frame dependence is 
minor and at present of no real theoretical concern; we 
do not document that finding, since the differences are 
only hardly seen in plots. 



3. Higher order contributions to the current operator in 
(/c/mjv) expansion 

In the standard calculational scheme the Siegert form 
of the current operator is used together with explicit 
meson-exchange contributions not accounted for by the 
Siegert part. The charge density operator in the Siegert 
part is of one-baryon nature and is taken to be nonrela- 
tivistic in the standard calculations. However, the one- 
baryon purely nucleonic charge density operator has rel- 
ativistic corrections of order (fc/mjv) 2 . Contributions to 
the nucleon-A transition charge de nsity and to the two- 
nucleon charge density, used in Ref. for calculation of 
trinucleon elastic charge form factors, are also included; 
both are of the relativistic order (fc/mjv) 2 . The resulting 
special relativistic corrections, taken into account in this 
paper, reduce the cross sections; they appear beneficial; a 
characteristic result is shown in Fig. El The effect shown 
there is dominated by the one-nucleon charge-density cor- 
rection; the two-nucleon charge-density shows noticeable 
effects in some spin observables, whereas the nucleon-A 
transition charge appears to be insignificant for all cal- 
culated observables of this paper. Correspondly large 
corrections of the same origin were also found in photo 
reactions on the deuteron Thus, the results of this 
subsection are not surprising. The current corrections of 
this subsection should be included in future calculations 
of e.m. reactions. 



V. SUMMARY AND CONCLUSIONS 

The paper improves our preliminary description of 
photo reactions in the three-nucleon system Q; the 
present description includes three-body photo disintegra- 
tion. The hadronic interaction is based on CD Bonn 
and its realistic coupled-channel extension CD Bonn + 
A. The initial and final hadronic states are calculated 
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FIG. 15: Differential cross section of proton-deuteron radia- 
tive capture at 200 MeV nucleon lab energy as function of the 
cm. nucleon-photon scattering angle. Results of the coupled- 
channel potential with A-isobar excitation derived from our 
standard approach (solid curve) are compared with results 
including relativistic one-nucleon charge corrections (dashed 
curve). The experimental data are from Ref. 25]. 



without separable expansion of the underlying interac- 
tion. The contributions to the e.m. current correspond 
to the hadronic interaction, though full current conser- 
vation could not be achieved. 



The paper isolates the A-isobar effects on the con- 
sidered observables. Besides A-isobar effects of effec- 
tive two-nucleon nature the A isobar yields an effective 
three-nucleon force of the Fujita-Miyazawa type and of 
the Illinois pion-ring type; meson exchanges other than 
pion exchange are included. The exchange currents me- 
diated by the A isobar are of effective two-nucleon and 
three-nucleon nature; they are structurally consistent 
with corresponding hadronic contributions; they are pre- 
dominantly due to the transition contributions of Fig. 
and due to the one-baryon part of Fig. 03 the diag- 
onal two-baryon contributions of Fig. connect small 
wave function components and are found to be quan- 
titatively entirely irrelevant. In the considered observ- 
ables the hadronic and the e.m. A-isobar effects are 
intertwined; they are not separated; their total effects 
are not very important given the scarcity of data, often 
still carrying large error bars. The A-isobar effects are 
more pronounced at higher energies; they are somehow 
smaller than the irreducible three-nucleon force effects of 
Refs. 0,0]; nevertheless, qualitatively both effects are 
quite similar. In contrast to the three-nucleon force ef- 
fects of Refs. fl^lT^ we see very small A-isobar effects at 
low energies; the reason, at least in part, is due to the dif- 
ferent choice of kinematics for correcting the theoretical 
failure in accounting for three-nucleon binding. 
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APPENDIX A: COUPLED-CHANNEL 
CURRENT OPERATORS 



Equation JQ| defines the general momentum space form 
of the e.m. current J M (Q) in the Jacobi coordinates 
of the three-particle basis. In contrast, this appendix 
gives its employed one-baryon and two-baryon parts, i.e., 
J"(Q) = J [11 ^(Q) + J [2lAt (Q), in respective one-particle 
and two-particle bases. We keep the three-momentum 
transfer Q and not the four-momentum transfer Q as 
independent variable since usually Qo is determined by 
the three- momenta of the involved baryons. Despite that 
strategy, Q 2 = Q 2 Q 2 , is taken to be zero in all e.m. form 
factors for photo reactions. The step from the single- 
particle representation of the current contributions to the 
three-particle Jacobi momenta is straightforward [2f| and 
not repeated here. The objective of this appendix is the 
definition of the used input for the current. 



1. One-baryon operators in nonrelativistic order 

The momentum-space matrix elements of the one- 
baryon current operator have the general form 



(k'B'|J W"(Q)|kB) = 6(k' - Q k)jW£(Q,k', k) 



pW w (Q,k',k) = e(Q 2 ), 
jW w (Q,k',k) = ^{e(Q 2 )[k' + k] 



(A2a) 



1 



J&(Q,kM0 2mA 

p W A (Q,k',k)= 5 I°(Q 2 
j AA (Q,k',k) 1 



;(Q 2 ) + K (Q 2 )][z<rxQ]}, (A2b) 
sS(Q 2 )[ 8 SxQ]T 2 , (A2c) 



(A2d) 



{ 3 I°(Q 2 )[k' + k] 



2mA 

+ 3a 1 (Q 2 )[^a x Q]}. (A2e) 

The nucleonic e.m. form factors of Eqs. l!A2jl are 
parametrized as linear combinations of the isospin- 
dependent Dirac and Pauli form factors fi(Q 2 ) and 
f2(Q 2 ), which at zero four-momentum transfer Q 2 are 
the charge and the anomalous magnetic moment of the 
nucleon, i.e., 

<Q 2 ) = ~ [A IS (Q 2 ) + A IV (Q 2 K] , (A3a) 

k(Q 2 ) = \ [ff{Q 2 ) + ti v (Q 2 )r z ] , (A3b) 

the superscripts IS and IV denote their isoscalar and 
isovector parts, respectively. The e.m. form factors re- 
lated to the A-isobar are parameterized [lffl l2ll by 



mA 



VAN 



(Al) 



m N (l + Q 2 /A 2 ANA ) 2 (l + Q 2 /A 2 AN ^/ 2 > 

(A4a) 

gl°(Q 2 ) = ^{fl s (Q 2 ) + fI v (Q 2 ) - Q 2 /(4m 2 r ) 

x[/f(0 2 ) + / 2 IV (Q 2 )]}i(l + r Az ), 

(A4b) 

The values of the parameters are /l«ajv = Ma = 

4.35 mat, [In being the nuclear magneton, Aajv,i = Aa = 
0.84 GeV and A AN 2 = 1-2 GeV. 



with k' (k) and Q being the final (initial) single-baryon 
momentum and the four-momentum transfer by the pho- 
ton, respectively, and B' (B) being TV or A depending on 
the baryonic content of the final (initial) state. All com- 
ponents of j^j^(Q,k',k) are still operators in spin and 
isospin space; the spin (isospin) operators of the nucleon, 
A isobar and the nucleon-A transition are denoted by a 
(t), <ta (ta) and S (T), respectively. The one-baryon 
charge density and spatial current operators, diagram- 
matically defined in Figs. HEl and used in the calcula- 
tions of this paper, are listed below: 



2. Two-baryon operators in nonrelativistic order 

The matrix elements of the two-baryon current opera- 
tor have the general form 

<k' 1 k' 2J B'|J,! 2lAl (Q)|kik 2 i?) 

= 5(k[ + k 2 - Q - ki - k 2 ) (A5) 

xi^ s (Q,k' 1 -k 1 ,k' 2 -k 2 ) 

with k- (k^ being the final (initial) single-baryon mo- 
menta; a denotes the exchanged meson or the two mesons 
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in case of nondiagonal currents); the baryonic contents 
B' and B being TV (A) correspond to the two-nucleon 

(nucleon-A-isobar) states. All components j'a^ B (Q, k^ — 
ki, k 2 — k 2 ) are still operators in spin and isospin space. 



The two-baryon spatial current operators, diagrammati- 
cally defined in Figs. G] - El and used in the calculations 
of this paper, are listed below: 



jiriwtQ.Pi.Pa) = - f™(Q 2 ){[iri x T 2 ] z F^ N (p 2 2 )(p 2 ■ 0-2)0-! + (1 «-> 2)} 

+ A IV (Q 2 )[iT! xx 2 ] z ^ s Ar (p?,p2)( Pl . (Tl )(p 2 . ( T 2 )(p 1 -p 2 ), (A6a) 

jJLrtQ.Pi.Pa) = - /r(Q 2 ){[iT! x T 2 } z F c p ° N n N (p 2 2 )[(cT 2 x pa ) x o-r] + (1 <- 2)} 

+ A IV (Q 2 )[iT! x r 2 ] 2 F^(p 2 ,p 2 )[( Pl x o-r) • (p 2 x o 2 )]( Pl - p a ) 
+ /i IV (Q 2 )[^i x r 2 ] z F^(p 2 , P 2 )( Pl - Pa ) 

- A IV (Q 2 )[^i x T 2 ] z F^ s N (p 2 l7 pl)Q x [(p a x o*i) x (p 2 x o- 2 )], (A6b) 
flwCQ.Pl.Pa) = - A IS (Q 2 ){(^i ■ x 2 )F^W(p 2 ,P^)(P2 • o- a )[i Pl x p 2 ] + (1 <- 2)}, (A6c) 
jL 2 Liv(Q ; Pi,P2) - - A IV (Q 2 ){t 22 ^ W (Pi,P2)(P2 • <r 2 )[ipi x P2 ] + (1 <- 2)}, (A6d) 



jl 2 L(Q,Pi,p 2 ) = - A IV (Q 2 ){[^i x T 2 ] z ^Sv(p|)( P2 • S 2 )o-a + (1 <- 2)} 

- A IV (Q 2 ){[*T! x T 2 } z F^ N (p 2 2 )(p 2 ■ o- 2 )S x + (1 <-> 2)} 

+ fl V (Q 2 ){[iTi x T 2 ] z ^X os / V (P 2 ) P2)(Pi -^i)(P2 • S 2 )( Pl -p 2 ) + 2)}, (A7a) 
jS JV (Q,Pi,p 2 ) = - A IV (Q 2 ){[^i x T 2 ] Z ^(P 2 )[(S 2 x P2 ) x ax] + (1 - 2)} 

- A IV (Q 2 ){NTi x x 2 ] z F p c ^(p^)[(<T2 x pa) x Si] + (1 ~ 2)} 

+ A IV (Q 2 ){^! x T 2 ] z i^-(p 2 ,p 2 )[( Pl x o-O • (p 2 x S 2 )](pi - P2 ) + (1 ~ 2)} 

- A IV (Q 2 ){[^i x T 2 ] Z F^°MP 2 ,P 2 )Q x [(pi x o-i) x ( P2 x S 2 )] + (1 «-> 2)}, (A7b) 
jS A iv(Q>Pi>P2) = - A IS (Q 2 ){(ti • T 2 )^ s AA r(p 2 ,p^( P2 • S a )[i Pl x pa] + (1 <-» 2)}, (A7c) 
jL 2 Liv(Q>Pi>P2) = - A IV (Q 2 ){7 1 2, J FL d ; s AJV (p 2 ,p2)(p 2 • S 2 )[ipi x p 2 ] + (!«-» 2)}, (A7d) 



j1 2 1a(Q.Pi,P2) = - fl V (Q 2 ){[irx x T A2 ] z F c ™ A d (p 2 2 )(p 2 ■ a A2 )a, + (1 «-> 2)} 

- A IV (Q 2 ){[^ai x r 2 ] z F™ n A d (p2)(p 2 • o-a)o- Aa + (1 «-» 2)} 

+ A IV (Q 2 ){^r x r A2 ] z ^ m A c A d (p 2 ,p2)( Pl • o-!)(pa • <x A2 )( Pl - P2 ) + (1 <-> 2)} 

- A IV (Q 2 ){[^ x T 2 ] 2 F- n A c (P 2 )(P2 • S 2 )S} + (1 ~ 2)} 

- A IV (Q 2 ){[^! x T 2 ],F^ A c (p 2 )(p 2 • Sl> + (1 ~ 2)} 

+ A IV (Q 2 ){[^1 x T 2 ],F^'(p?,Pa)(Pi • S+Xpa • S 2 )( Pl - p 2 ) + (1 - 2)}. (A8) 



We note: The contribution to the two-nucleon 
p-exchange current, proportional to F„jyjy (pf ,Pa) 
Eq. l|A6bll is not contained in the standard collection of 
exchange currents of Refs. [fil H3, Il7l l2fj. used by us till 
now in the context of other potentials; it is necessitated 
in this paper by the full form of the p exchange imple- 
mented in the CD-Bonn potential. Other contributions 
arising from the full p exchange jj3] are of higher order 



compared to -FJjvat (p 2 , p 2 ) and therefore are neglected 
in our standard calculations; their effect is discussed in 
Sec.lTvTm 

The -F-functions used in the above expressions are 
potential-dependent. For meson-exchange potentials 
they are built from meson-baryon coupling constants, 
hadronic form factors and meson propagators. For con- 
tact currents the F-functions have the following forms: 
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K°nn(p 2 ) 



&-K 2 m? N 4ir m 2 + p 2 ' 



(A9a) 



peon ,„2n _ 1 gp(l + /p/gp) 2 ^pW(P 2 ) 

pN N \P ) ~ 87r 2 TO 2^ 47r m 2 +p 2 • 



F C P 1\(P 2 ) 



peon, d/ 2\ 



(A9b) 

1 9^ fnNA ^ 7 ttNn(P 2 )^An(P 2 ) 



8ir 2 m z N 4?r f nNN m 2 + p 2 

i g 2 (i + /p/g P ) 2 f P NA 

8ir 2 m 2 N Air f pNN 

F p nn{p 2 )F p an(p 2 ) 



(A9c) 



(A9d) 

1 gl /ttAA ^ r 7rWAr(p 2 ).? r 7rAA(p 2 ) 



m 2 + p 2 



8ir 2 m 2 N 4?r /^atjv 



m 2 + p 2 



(A9e) 

r-icon, c/ 2\ J- g 2 f-wNA •^rAjvjP) Mnfi 

nAA (P ) - ^T^/ 2 ^ m 2 + p 2 • 

For meson in flight currents the corresponding expres- 
sions are 



1 



^STuCPi. P2) = - -5^-3 [Kb'b(p 2 i) - Ks'siPl)] , 

Pi P2 

(AlOa) 



TTimcsl,' 2 2\ 

r p NN\Pl,P2) - 



Am 2 N 



(1 + /p/5p) 2 PNN 



pmes ( 2 2\ 
fpNN{Pl>P2), 



pmes, d(e)/ 2 „2\ 1 
^ttAA lPl)P2i - ~ 2 _ „2 
Pi P2 



x[CA d(e) (P?)-CA d(C) (P^)]- 

(AlOc) 

Finally, the functions for nondiagonal meson-exchange 
currents (also called dispersion currents) are defined to 
be 



F dis / 2 2n _ i g a gp m N 

KpNN{Pl,P2) - in2m 2 N in — 



9a/3j 



FgNN (Pi) FpNNjPz) 

m 2 a + pi mj| + p\ 

1 gaff/3 /ttJVA "T-tv 



, (Alia) 



-Fa/3AAr(Pli P2) — T~2 T~ i 7 #"07 

47T 2 m^ 47T /jrJVAT ™a 

^aJVJv(Pl) •? 7 /3AA r (P2) 



m 2 + p\ ml + p\ 



(Allb) 



The meson-nucleon coupling constants <jr Q and f p are 
listed in Table I of Ref. whereas other hadronic pa- 
rameters, i.e., coupling constants f a B'B, meson masses 
m a and hadronic form factors F a B'B{p 2 ) are those of 
Ref. The e.m. meson-photon coupling constants 

have the values g p7r7 = 0.56 and g W7r7 = 0.68 according 
to Ref. EH- 



3. Operator corrections of lowest relativistic order 

Sample operator corrections of relativistic order for the 
charge density are given. They are of one-baryon and of 
two-baryon nature: 



N 



{Q 2 + [i<rx (k' + k)] -Q}, 



(Al2a) 



p^ rc (Q,k',k)= - 



1 



-9an(Q 2 ) 



[2] rc v 
PUw(Q,Pl,P2) 



4mNlTlA' 

x [iS x (k' + k)] QT Z , (A12b) 
^[/ 1 IS (Q 2 )r 1 -T2 + / 1 IV (Q 2 )r 2 ,] 

X i&(p!)(<>-l • Q) (0-2-P2) 



+ (1 <-> 2). (A12c) 

The contributions llA12|l are the Darwin-Foldy and spin- 
orbit corrections of the one-nucleon charge density, the 
one-baryon correction due to nucleon-A transition and 
the two-nucleon correction due to it exchange, respec- 
tively; the two-nucleon contribution l|A12cjl is local 
and therefore often exclusively used; there are however 
other nonlocal two-nucleon contributions of the same or- 
der. The contributions l|A12j) are used in this paper in 
Sec. II V C3l for the Siegert form of the current. Since they 
are relativistic corrections, they violate current conser- 
vation in the considered order. However, the calculated 
trinucleon elastic charge form factors need all three con- 
(AlOb) tributions in order to become almost quantitatively con- 



sistent with the experimental data 



3t qt 

m 



APPENDIX B: INTEGRAL EQUATION FOR 
CURRENT MATRIX ELEMENT 

This appendix calculates the current ma- 
trix elements of two- and three-body photo 
disintegration of the trinucleon bound state, 



i.e., (V4 _) (q/) 

i/ Q/ |j A '(k 7 ,K + )e AI (k 7 A)|B) 

(V^ ) ( P/ q / )^o / |jnk 7 ,K + )e A1 (k 7 A)|i3). 

The antisymmetrized fully correlated three-nucleon 
scattering states of internal motion in nucleon-deuteron 
channels, i.e., (ipa \ l '\}) v a s \, and in three-body breakup 

channels, i.e., (ipQ \pf<lf)vof\, are not calculated explic- 
itly; they are calculated only implicitly when forming 
current matrix elements. We introduce the state \X(Z)), 
defined according to 



and 



\X{Z)) = (l + P)^(k 7! K + )e M (k T A)|B) 
+ PT(Z)G (Z)\X(Z)), 

oo 

\X(Z)) = Y / [PT(Z)G (Z)r 



(Bla) 



n=0 



(l + P)j' i (k 7 ,K + )e /i (k 7 A)| J B), (Bib) 
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as intermediate quantity with Z = Ei + iO being the 
three-particle available energy and T{Z) being the two- 
baryon transition matrix. Equation ijBlal is an integral 
equation for \X(Z)), analogous to that for the multichan- 
nel transition matrix U(Z) of Ref. jg: Both equations 
have the same kernel, only their driving terms are differ- 
ent. We therefore solve Eq. IIBlaJl according to the tech- 
nique of Ref. [3| , summing the Neumann series fBlbj) for 
\X(Z)) by the Pade method. Once \X{Z)) is calculated, 
the current matrix elements required for the description 
of two- and three-body photo disintegration of the trin- 
ucleon bound state are obtained according to 



^ ) (q/)^ / l^(k 7 ,K + )e Al (k 7 A)|i?) 
1 -.(M<lf)v af \X(Z)), 



(B2a) 



(^- ) (p / q / )^ 0/ |^(k 7 ,K + ) ep (k 7 A)|i?) 

= -^<0 o (p/q/)^o / |(l + P)[/ I (k 7 ,K + ) eA1 (k 7 A)|i?) 



T(Z)G (Z)\X(Z)}}. 



(B2b) 



The current matrix element required for the description 
of radiative nucleon-deuteron capture is related to that 
of two-body photo disintegration by time reversal as de- 
scribed in Ref. §j. When calculating total two- and 
three-body photo disintegration cross section a, the inte- 
gration over all final states can be performed implicitly, 
i.e., 



x 6(Ei -H - JJ J ).7' i (k 7 ,K + )e M (k 7 A)|B}, (B3a) 



7 M B \ 



£ Im{( J B|^(k 7 ,K + )e /1 (k 7 A)]t 

MbA 

G(^+ i 0)/ l (k 7 ,K+)e A1 (k 7 A)| J B)}. (B3b) 



The auxil iary state G{E l + iO)^(k 7 , K + )e M (k 7 A)|B) of 
Eq. l|B3hl) is related to \X(Ei+iO)) according to 



G^+iO^k^K+^^A)^) 

= i (1 + P)G {Ei +*0) [f (k 7 , K+)e M (k, A) \B) 
+ T(E i +iO)G {E i +iO)\X(E i +iO))] . (B3c) 



The total lab cross section is then obtained in the form 



°= E I m {(S|[^(k 7 ,K + ) £ ,(k 7 A)]t 

x (l + P)G (E i +iO)[j"(l h ,K + )e fi (l h X)\B) 

+ T(E i +iO)G {E i +iO)\X(E i +iQ))]}. (B3d) 



We note that Eqs. l|B3jl performs the integration over all 
final states implicitly using the nonrelativistic Hamilto- 
nian in contrast to the strategy of Eq. ijlO)) . 
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